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Excitation of slow MJO-like Kelvin waves in the equatorial 
atmosphere by Yanai wave-group via WISHE-induced 
convection 
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The intraseasonal Madden- Julian oscillation (MJO) in- 
volves a slow eastward-propagating signal in the tropical 
atmosphere which significantly influences climate yet is not 
well understood despite significant theoretical and observa- 
tional progress. 

We study the atmosphere's response to nonlinear "Wind 
Induced Surface Heat Exchange" (WISHE) forcing in the 
tropics using a simple shallow water atmospheric model. 
The model produces an interestingly rich interannual behav- 
ior including a slow, eastward propagating equatorial west- 
erly multiscale signal, not consistent with any free linear 
waves, and with MJO-like characteristics. It is shown that 
the slow signal is due to a Kelvin wave forced by WISHE 
due to the meridional wind induced by a Yanai wave group. 
The forced Kelvin wave has a velocity similar to the group 
velocity of the Yanai waves, allowing the two to interact 
nonlinearly via the WISHE term while slowly propagating 
eastward. These results may have implications for observed 
tropical WISHE-related atmospheric intraseasonal phenom- 
ena. 



1. Introduction 

The "Madden Julian Oscillation" (MJO) [Madden and 
Julian, 1971, 1994; Zhang, 2005] is the most significant in- 
traseasonal phenomena in the tropical atmosphere. It cou- 
ples convection and large-scale circulation, and propagates 
eastward at about 5 ms _1 , through the Indian and East to 
Central Pacific oceans. This is a surprisingly slow propaga- 
tion when compared with a typical moist Kelvin wave speed 
(15-20 ms- 1 ). 

One of the mechanisms suggested to explain the MJO was 
wind induced surface heat exchange (WISHE) [Emanuel, 
1987; Neelin et al, 1987], according to which air-sea latent 
heat fluxes enhanced by surface winds are of paramount 
importance for forcing and organizing tropical convection. 
Emanuel [1987] and Neelin et al. [1987] found Kelvin-like so- 
lutions due to the convection-wind coupling in linear WISHE 
models. While these early solutions seem to differ signifi- 
cantly from observed MJO characteristics [e.g. Zhang, 2005], 
recent observational and model evidence suggests an impor- 
tant role for WISHE in MJO dynamics [Moloney and So- 
bel, 2004; Sobel et al., 2008]. In addition, nonlinear WISHE 
effects seem to be important in simulated MJO-like distur- 
bances [Raymond, 2001; Moloney and Sobel, 2004]. 
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In light of these findings we investigate an equivalent- 
barotropic shallow water atmospheric model with a simple 
nonlinear WISHE formulation, on a beta plane. The nu- 
merical experiments yield intriguing dynamics, including a 
slow eastward propagating signal that shares some of the ob- 
served MJO characteristics, and its mechanism is analyzed 
below and also presented in the thesis of Solodoch [2010] 
(hereafter S2010). Given the simplicity of the model, it is 
not possible to conclusively demonstrate that this mecha- 
nism is applicable to the observed MJO, but we feel it is 
nevertheless a most interesting and nontrivial mechanism 
which can enrich our understanding of tropical dynamics. 

2. The model 

The model equations are the shallow water equations with 
a WISHE heating forcing term (Q) [Gill, 1980], 
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where u and v are the zonal and meridional perturbation ve- 
locities, rj is the perturbation (hydrostatic) pressure, v the 
eddy diffusivity coefficient, and r a Newtonian cooling time 
scale. In the standard experiment, mean surface easterlies, 
U = Ux, which tend to prevail in much of the equatorial 
Pacific are added to the model velocity in the WISHE term 
only, with U = — 2 m/s. We also set r = 5 days, rj aa t = —5 
mbar and A = 2 • 10~ 6 m -1 . The equivalent depth is set 
to H = 50 m, consistent with observed convectively cou- 
pled waves in the tropics having H = 10 — 60 m [Wheeler 
and Kiladis, 1999]. The typical tropical balance between 
convective heating and adiabatic cooling is thus implicit in 
(l)-(3) through this choice of H, and Q represents only the 
positive-definite convective heating caused by surface heat 
flux variations. The model is periodic in the east-west coor- 
dinate (x), extends from 60S to 60N and 25,600 km in the 
east-west direction, and is based on the MIT GCM [Marshall 
et al., 1997a, b]. 



3. Results: slow eastward propagation of 
equatorial westerlies 

A Hovmoller diagram of the zonal velocity u at the equa- 
tor, for days 80-160 after the initialization, is shown in 
Fig. la. The most prominent feature is an eastward propa- 
gating disturbance, significantly slower than the equatorial 
Kelvin wave (22 m/s in this model). The slow signal is an 
irregular, a wide band, multiscale phenomenon in space and 
time, and the propagation speed varies between 8-12 m/s. 
One wonders, of course, if this slow eastward propagating 
equatorial disturbance is of any relevance to the MJO, but 
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Figure 1. Hovmoller diagrams at the equator from exp 
1. Top: zonal velocity, u [m/s], with lines showing propa- 
gation speeds of 22 m/s and 9.5 m/s. Middle: meridional 
velocity, v [m/s]. Bottom: Spectral power of u(x,y,t), in 
log scale, calculated by Fourier transforming in (a;, t) at 
each latitude, and then summing the power between ±30 
degrees latitude. Zonal wave number is for an equatorial 
circumference of 25,600 km. 

understanding its dynamics seems a worthwhile goal in any 
case. 

Interestingly, the anomalous zonal wind is not symmetric, 
and a histogram analysis (not shown) shows that easterly 
anomalies are usually weaker than the westerly anomalies. 
This is reminiscent of "Westerly wind bursts" observations 
over the tropical oceans [Harrison and Vecchi, 1997] associ- 
ated in some cases with MJO events. 
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Figure 2. The meridional structure of the slowly prop- 
agating signal. The blue line shows a meridional cross 
section showing the root time-mean square of zonal ve- 
locity anomalies at the middle meridian of the domain. 
Green: a Gaussian meridional structure of the Kelvin or 
v velocity of the Yanai wave. Red: a Gaussian fit based 
on a Rossby radius of deformation smaller by a factor of 
\/2 than the model value. 

A Hovmoller diagram of the equatorial meridional veloc- 
ity, v, (Fig. lb) shows an eastward propagating wave group, 
propagating with the same speed as the slow signal seen in 
the u Hovmoller diagram. The v wave group envelope is 
highly coherent with the slowly eastward propagating u sig- 
nal, while showing complex patterns of phase propagation 
within the eastward propagating envelope. The v Hovmoller 
and the spectrum shown in Fig. lc show that the v signal 
corresponds to a group of Yanai waves propagating with the 
slow signal. It will be shown below that the speed of propa- 
gation of the slow signal is indeed compatible with the group 
velocity of long Yanai waves, and the interaction between the 
u and v fields will be explained. 

Dispersion curves of equatorial waves are superimposed 
over the contours of the power spectrum in Fig. lc, and most 
of the spectral power is seen to be concentrated near these 
dispersion curves. The linear wave modes with the most 
power are the Yanai wave and the Inertia-Gravity n = 1 
wave, and the Kelvin wave is also excited. 

However, significant power is concentrated in a feature 
not associated with any linear wave mode, which is located 
in k — oj space beneath the dispersion curve of the Kelvin 
wave. This feature corresponds to the slowly propagating 
signal seen in the Hovmoller diagrams. A cross section along 
the slowly propagating signal (Fig. 2, blue line) shows its 
variance to be symmetric in latitude with a Gaussian struc- 
ture centered around the equator, similarly to a Kelvin wave. 
It also lies close to a linear line in k — ui space, and appears 
in the t) field, yet is virtually non-existent in the v power 
spectrum (S2010). These are all characteristics of a Kelvin 
wave, except that this feature has the "wrong" meridional 
scale when compared to the Gaussian Kelvin structure given 
the model's equivalent depth (Fig. 2) and this is further ex- 
plained below. We also note that most of the energy of the 
"sub-Kelvin" signal is in the k — 1 planetary scale, and the 
energy decreases with k. The sub-Kelvin variability explains 
50% of the zonal wind variability at the equator. The above 
description makes it clear that the slow eastward propagat- 
ing signal clearly combines both low frequency Kelvin-like 
waves, and higher frequency Yanai waves. 

4. The mechanism 
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Figure 3. Demonstrating the longitudinal phase between the heating and the slow signal wind response. Arrows: winds 
corresponding to the slowly propagating signal (obtained by filtering the full solution in k, ui space, keeping only the area 
under the kelvin wave dispersion curve corresponding to the show signal alone. Contours: the 2-day running time-averaged 
Yanai wave contribution to the heating {v Yanai }. 



We now show that the slow eastward propagating signal 
is a forced Kelvin wave, excited by the envelope of a Yanai 
wave-group. The efficient excitation is possible due to two 
factors. First, the forced Kelvin velocity is very similar to 
the Yanai group velocity. Second, the meridional structure 
of the slow Kelvin waves signal is identical to that of the 
square of the Yanai wave, which allows the Kelvin wave to 
be excited through the nonlinear WISHE heating term. 

Starting from the heating forcing term, note that the 
power series expansion of the WISHE term Q is, 



(ri+\ri*\)V(u + U) 2 +v 2 = \ri*\\U\ + V \U\ - 
k*l 2 1 2 , IV 



1 2 . VI | 2 . 



Although the formal convergence criterion for this power 
series is not always satisfied in the model, a correlation 
analysis of the different terms in this expansion with the 
full heating Q shows that the v 2 term accounts for much 
of the heating variability, v itself is, in turn, dominated 
by the contribution due to the Yanai waves. Now, we 
are interested in the slow forcing of a Kelvin wave by the 
meridional velocity due to the Yanai waves, vy ana i, lead- 
ing us to examine the time average of (v Yanai ) 2 . This 
represents the energy of the Yanai waves which therefore 
moves at the Yanai group velocity c g (k) and may be written 



(VYanai) = Efc H i V y) 



■J 2 Y ) k e ik ^- 



,t) 



+ c.c. 



where 



complex conjugate, and (v Y )k is the amplitude of the enve- 
lope of v 2 due to the Yanai waves at wavenumber k. 

An important insight here is that the Yanai v field has 
a Gaussian structure, = exp(— y 2 / {2R 2 ,)), so that its 
square results in WISHE forcing of the form Q~u 2 ~ 
exp(-2j/7(2i^)) = exp(-y 2 /(2Ri)), where R x = R /V2. 
This also implies an equivalent depth of the forced signal 
of H/4 where H is the model's specified equivalent depth. 
In other words, the square of the Yanai wave v velocity has 
the Gaussian structure corresponding to a smaller equivalent 
depth, forcing a Kelvin wave like response with a similarly 
smaller effective equivalent depth. This implies that this 
forced Kelvin wave propagates slower than the free Kelvin 
waves in the system, leading to the observed slowly propa- 
gating signal. 

We found in the numerical simulation above that Q ex- 
cites a Kelvin-like signal with a very small v velocity. Con- 
sider therefore the following set of equations for a forced 
Kelvin wave by setting v — and ignoring dissipation, 



u t = 
/3yu = 



Tit + Hu x 



-Vv, 
-G(«? 
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ik(x- 

2R 2 e K 
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(5) 
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Assuming G{v Y )k to be a constant amplitude of the forcing 
at wavenumber k (see discussion below), the forced Kelvin 
solution would be a superposition of the solutions with the 
same functional form as the forcing, 



ik(x- 



Cgt) 



(m ,?7o) 



(8) 



Inserting (8) into (5-6) leads to the condition c g = /3R 2 — 
Co/2 where Co is the free Kelvin wave speed in the model. 
Interestingly, this condition is satisfied exactly for a Yanai 
wave of fc = as can be verified from the Yanai disper- 
sion relation, and more generally is approximately satisfied 
for long Yanai waves of — 5<k < 10. This interesting co- 
incidence allows the slowly propagating solution (8) to be 
excited by the Yanai waves. Inserting this relation into (7) 
yields the forced wave solution, 



u 



iG{v 2 Y ) k 
/c(cq — c\ 



r)o = c g u , 



(9) 



These expressions are found to be consistent to a very good 
approximation with the numerical simulation. The dis- 
turbance clearly has the strongest response at the longest 
wavenumbers (k = 1), as is actually found for the sub-Kelvin 
spectrum in our simulation (and, interestingly, also for the 
observed MJO). Fig. 2 shows a meridional cross section 
through the slowly propagating signal (blue), demonstrat- 
ing that its structure is indeed very nearly a Gaussian (red) 
corresponding to the square of the Yanai wave structure 
(green). This provides an explicit and strong evidence of 
the nonlinear forcing mechanism described above. 

Observations show that the MJO convection is mostly 
co-located with its surface westerlies in the west-Pacific and 
with surface convergence in the Indian ocean [Zhang, 2005]. 
The Hovmoller plot of u and v (Fig. la,b) show that the 
maximum of the Yanai envelope (and thereby its contribu- 
tion to WISHE) is nearly co-located with the maximum of 
the sub-Kelvin signal at all times. Similarly, the snapshot 
of Fig. 3 shows the maximum Yanai induced v 2 WISHE en- 
velope to be co-located with, or slightly lead, the maximum 
sub-Kelvin westerlies, similar to the corresponding relation- 
ship between MJO convection and circulation, as observed. 

The Yanai waves which force the slow signal may them- 
selves be destabilized by the WISHE term in the presence of 
easterly winds [Emanuel, 1993]. However, the Yanai fields 
on the left hand side of the continuity equation (3) (77, u x , 

v y ) are antisymmetric in y (e~2( H//fl ' 2y/R) and therefore 
cannot be forced by the symmetric, Gaussian, Kelvin wave 
fields (u and rf) even when raised to some power, justifying 

the above approximation of G{v Y )k as a constant amplitude 
of the forcing not affected by the Kelvin waves. 
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5. Sensitivity study 

The main features of the slowly propagating signal are 
found to be robust with respect to different meridional 
boundary condition, initial conditions, and the value of the 
mean parametrized zonal velocity, U (as long as it is east- 
erly). The signal also appears if the WISHE forcing term is 
limited to the equatorial domain only, as long as this forcing 
domain is not smaller than about 6 degrees latitude. The 
slow signal is also quite robust to the values of model param- 
eters, yet does not appear, of course, for large values of the 
dissipation parameters or too small values of the WISHE co- 
efficient. Given that the WISHE coefficient, A, idealizes the 
product of a transfer coefficient and several other poorly 
constrained thermodynamic factors in a bulk transfer for- 
mula [Emanuel, 1987; Neelin et al., 1987, e.g.,], it is not 
obvious whether the value used here is realistic. Additional 
sensitivity experiments and analysis is given in S2010. 

6. Conclusions 

An equivalent-barotropic shallow water model driven by 
nonlinear WISHE convective heating was used to study trop- 
ical atmosphere intraseasonal phenomena. A slowly east- 
ward propagating signal appears robustly in the experi- 
ments, with about half the speed of a linear Kelvin wave. 
The slowly eastward propagating equatorial signal is inter- 
esting because it shares some major characteristics with the 
large scale fields of the MJO, some of which are not ex- 
plained by current theories. 

The slow eastward propagating signal was shown to be 
a forced Kelvin wave with a meridional spatial structure 
consistent with an equivalent depth smaller than that spec- 
ified in the model. The signal is also coherent with a group 
of long Yanai waves which has a group velocity similar to 
that of the slow forced Kelvin wave. We showed that the 
Yanai group nonlinearly forces the slow Kelvin-like signal 
through the WISHE term, and explained how the nonlinear 
WISHE term leads to the smaller effective equivalent depth 
and slow propagation of the signal. The suggested mecha- 
nism predicts the slower Kelvin wave meridional structure 
and speed, as well as the right order of magnitude of the am- 
plitude. The excited signal in the numerical simulations is 
strongest at k = 1 and 2, again consistent with the suggested 
mechanism. The phenomenon has a multiscale character in 
space and time, and an interesting and novel aspect of the 
mechanism proposed here is the critical role of higher fre- 
quency motions involving the Yanai waves. Concurrently 
with Solodoch [2010] and the present study, a somewhat 
related idea has recently been presented in a conference ab- 
stract [Yang and Ingersoll, 2009]. 

The slow signal is reminiscent of some of the main fea- 
tures of the MJO, in particular its wide-band nature, slow 
eastward propagation, its Kelvin-like large scale structure 
at and east of its convection center, and its spectral power 
peaks at the k — 1 and 2 wavenumbers. Additionally, there 
is asymmetry in zonal wind toward higher westerlies than 
easterlies, particularly close to the convection center. 

Preliminary results to be reported elsewhere [Solodoch 
et al, 2010] show that the mechanism found here is also 
relevant to more detailed models, with significantly more re- 
alistic convective parameterizations. Although it is clearly 
premature to propose that the mechanism found here is ap- 
plicable to the MJO, exploring the suggested mechanism in 
more realistic models and in observations may prove fruitful. 
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